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Tree models for rigidity percolation are introduced and solved. A probability vector describes 
the propagation of rigidity outward from a rigid border. All components of this "vector order 
• parameter" are singular at the same rigidity threshold, p c . The infinite-cluster probability Poo is 

usually first-order at p c , but often behaves as Poo ~ AP 00 + (p—pc) 1 ^ 2 , indicating critical fluctuations 
superimposed on a first order jump. Our tree models for rigidity are in qualitative disagreement 
with "contraint counting" mean field theories. In an important sub-class of tree models "Bootstrap" 
percolation and rigidity percolation are equivalent. 
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I. INTRODUCTION 



Soon after the resurgence of interest in percolation phenomena, the elastic constants of depleted materials were 
studied. Although early work suggested! that the conductivity and elasticity exponents were the same, it was soon 
realised that the elasticity exponents were usually differenfn and in particular one must draw a distinction between the 
elasticity of systems which have only "central forces"H and those which also have "bond-bending" forces. If a system 
has bond-bending forces, the percolation geometry is in many way similar to that of the connectivity percolation 
problem. Of interest in this paper is the fact that when a system is supported by only central forces, the percolation 
' geometry is very different than that occuring in connectivity percolation. We illustrate this difference by developing 
t-H , and solving models for rigidity percolation on trees and by comparing those models with the analogous results for 
connectivity percolation on treesQ. Many of the concepts we develop using tree models can be extended to regular 
lattices, as will be elaborated upon in the paper. 

There have been several different groups of scientists and engineers interested in the ability of central force structures 
to transmit stress. Besides its intrinsic interest, this problem is relevant to the analysis of engineering structures, 
glasses, granular materials and gelsQ. The straightforward way to study this problem is to construct particular models 
which have only central forces and to study the types of structures which support stress. In the physics community, the 
standard model has been lattices composed of Hooke's springs. Direct solution of the force equations for these lattices 
has provided quite variable estimates of the percolation threshold, and considerable controversy about the critical 
exponentsETlj. In the mathematics commupitjVithere has been a long history of attempts to related the-connectivity 
of a "graph" to its ability to support stresst3~ll3. The majority of physicists were unaware, till recentlyOllil, that there 
is a rigorous theorem which relates connectivity to rigidity but only for planar graphs. Of more practical importance 
is the fact that there is a fast algorithm^ by which this theorem can be used to actually find the infinite clusterta 
and stressed backboned of planar graphs (e.g. the triangular lattice with central forces). These results are relevant to 
random latticestJ, which are in many cases of most practical interest. 

There are two different types of mean field theory available for the rigidity transition. The first, hased on an 
approximate "constraint counting", predicts a second order transition in the "number of floppy modes" t£l, and has 
been extensively applied to the rigidity of glasses and gels. However it was realised in that paper and recently 
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quantitatively confirmedE-3 that the number of-floppy modes per site does not appoach zero at the percolation point. 
There is also a recent continuum field theorytJ which predicts first order rigidity, but the connection between the 
model parameters and the lattice parameters relevant to rigidity are not clear in that analysis. The tree models 
developed here provide a more complete mean field theory for the rigidity transition. We find that the number of 
floppy modes is continuous near the rigidity transition, but that the infinite cluster probability shows a first order 
jump. We also find that one sub-class of our tree models are equivalent to tree models for bootstrap percolationEZI, 
although on regular lattices they are not equivalent. 

The paper is arranged as follows. In the next section (Section II), we introduce the tree geometry and the vector 
probability (order parameter) used to describe the transmission of rigidity from a rigid border. Section III contains 
the detailed analysis of the tree models for both site and bond dilution. In Section IV we discuss the mechanism for 
first order rigidity and discuss the failings of the traditional constraint counting mean field theories in the light of the 
tree results. We also calculate the number of floppy modes, and show that even the second derivative is non-singular 
on trees. This is not too surprising, since surface bonds dominate on trees. Section V contains a brief summary and 
conclusion. 



II. THE GEOMETRY AND DEFINITION OF VARIABLES 



The structure of the tree models we consider is illustrated in Fig. 1. Following normal convention, we define z to 
be the number of branches of the tree (for example in Fig. la z = 5). In Fig. la, each site of the tree is connected 
by only one bond to a neighbouring site. In general we may have b bonds connecting neighbouring sites (for example 
in Fig. lc, 6 = 3). Thus two variables in our analysis arc z and b. A third important variable, g, is the number of 
degrees of freedom per site and is discussed in the next paragraph. The feature of the tree geometry which makes the 
analysis tractable is that we can calculate the probability of rigidity along separate branches of trees independently, 
and then combine the branches of the tree to form the final Bethe lattice. For example one "branch" of the tree of 
Fig. la is presented in Fig. lb. We use the letter P, with various subscripts, for the site probabilities of the entire 
tree (e.g. Fig. la), while we use T, with various subscripts, to denote the site probabilities of the branches of the 
trees. The qualitative behavior of T and P are the same, and we concentrate for the most part on the analysis of T. 

Each node (the sites in Fig. la,b and the ellipses in fig. lc) represents a "joint" (a point-like node) or "body" (see 
below) on a lattice or "graph" , and is assigned a certain number of "degrees of freedom" . In connectivity percolation 
each node is either connected or disconnected, so it has only one possible "degree of freedom", i.e. if a site is 
disconnected it has one degree of freedom, while if it is connected it has no degrees of freedom. If we consider a 
lattice of joints connected by central force springs, then each free joint has two translational degrees of freedom in 
two dimensions and three degrees of freedom in three dimensions. However when we make rigid clusters, they are 
rigid "bodies" so they also have rotational degrees of freedom. For example, a body in two dimensions has 3 degrees 
of freedom (two translations and one rotation), while a body in three dimensions has 6 degrees of freedom (three 
translations and three rotations). In general, we allow each site to have g degrees of freedom. Some practically 
important values for g are as follows, 

g = 1 for connectivity percolation, (la) 

g = d for a joint, (lb) 

g = d(d + 1) /2 for a body. (lc) 

Here d is the spatial dimension. We consider growing clusters from a rigid boundary drawn around the outer perimeter 
of the tree. If a rigid cluster grown from this boundary continues to grow indefinitely, we are above the percolation 
threshold, if it dies out we are below the percolation threshold. The behavior on crossing the percolation threshold 
depends on whether the transition is first order or second order, as will be discussed further below. In the case 
of connectivity percolation, there is only one degree of freedom per node, and we only have to keep track of the 
probability that connectivity is transmitted away from the boundary. In the case of rigidity percolation it is necessary 
to consider a larger set of site probabilities. In fact, each site may have 0, 1, 2, ...g degrees of freedom with respect to 
the boundary, so we define the probabilities Pq 7 .--Pg to be the probabilities that a site has between and g degrees of 
freedom (DOF) with respect to the boundary (a similar definition applies to the branch probabilities T). For example 
if 9 = 3 

(1) P?,{or T3) is the probability that a node has 3 DOF w.r.t. the border 

(2) Pi(or T2) is the probability that a node has 2 DOF w.r.t. the border 

(3) P\(or Ti) is the probability that a node has 1 DOF w.r.t. the border 

(4) Po(or Tq) is the probability that a node has DOF w.r.t. the border. 
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The vectors P and T act as order parameters for the rigidity percolation problem on trees. However, it is also 
possible to define these quantities on regular lattices and it is likely that an algorithm could be developed based on 
these probabilities. In fact for the case of a "diode response", a transfer matrix could be used - this would be a 
"directed rigidity percolation" and might be appropriate for granular media, where contacts only support compressive 
forces. 

In many physical problems, it is important to disinguish between a site which is overconstrained or stressed, Pb, 
and one which is rigid but not stressed (which has probability Pd). In particular, we have previously defined'' 
Poo — Pq = Pd + Pb toJae the "infinite rigid cluster" probability. This is closely analogous to the infinite cluster in 
connectivity percolationtS. In this analogy, the overcontrained or "stressed" bonds are analogous to the "backbone" 
in connectivity percolation. Also, just as the dangling ends in connectivity percolation carry no current, the dangling 
ends in rigidity percolation carry no stress. However for trees we found it clearer to first concentrate on Pinfmity-, so 
in this paper we do not discuss Pb ■ 

III. DILUTED BETHE LATTICES 

Consider Bethe lattices of co-ordination number z as shown in Fig 1. In general our parameters are g (the number 
of degrees of freedom per node), z (the co-ordination number - Actually we shall usually use a = z — 1), b (the 
number of bonds connecting each pair of nodes) and p (the probability that a site or bond is present). We first do the 
calculations for a branch of the trees (see Fig. lb. for a b = 1 case) and then join the branches together. To illustrate 
the method we first do the case b — 1, as illustrated in Figs. la,b with site dilution. 

A. Site diluted Bethe lattices with b=l 

On any tree, rigidity can only be transmitted to higher levels of the tree if there are enough bars present to offset 
the number of degrees of freedom of a newly added node. For connectivity percolation only one bar is needed. If a 
node is added to a g=2 tree, two bars are needed to offset the two degrees of freedom of the added node. In general, 
if a node with g degrees of freedom is added, rigidity is transmitted to the next level of the tree provided the node is 
occupied and provided at least g of the lower level nodes to which the added node is connected are rigid. We define 
the probability that a node is rigid to be To- The branch probabilities Tk with k — 0, l...g are then given by, 




s-i 

T g = l-J2?i (2) 

1=0 

The left hand side of Eqs. (2) refer to a node at the one higher level than the nodes on the right hand side. Since we 
are looking for asymptotic probabilities a long way from the rigid boundary, we expect the probabilites T/ to approach 
steady state values upon iteration of Eqs. (2). Similar expressions to Eqs. (2) are found when the transition is made 
from the branch probabilies T; (see Fig. lc) to the tree probabilities Pi (see Fig. la), except that we now combine z 
branches instead of z — 1 branches. Thus we find, for example 

Po=pY J {^jTl{l-T y- 1 (3) 
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In fact once we have solved the first of Eq. (2) and have found To, all of the other components of P and T follow. In 
particular, if To is first order at a particular p c , then all of the other components of T and P are first order at the 
same p c . Thus we concentrate on the behavior of T). ,— . 

It is interesting to note that Eq. (3) is the same as Eq. (2) ofhH which treats bootstrap percolation on trees (with the 
change of variables R = 1 — P, and g — m and / = z — to). In bootstrap percolation one considers that ferromagnetic 
order is propagated only if each site has at least to ferromagnetic neighbours. If we start with a ferromagnetic border, 
it is clear that Eq. (3), with the above change of variables, describes the propagation of ferromagnetic order outward 
from the border. The correspondence between bootstrap percolation and rigidity percolation is not exact on regular 
lattices, and it is not clear how to distinguish between these two cases in a continuum field theory calculation. 

Now we do some detailed solutions to the Eqs. (2). First we treat some simple solvable cases. 

Connectivity percolation (g = 1) n 
In this case the first of Eqs. (2) reduces to that found previouslyQ. For example for a = 3 

To = p(3T (l - T ) 2 + 3T 2 (1 - T ) + T 3 ) (4) 

which yields the trivial solution T = 0, and the non-trivial solution 

3- V(4/p-3) 

To = 2 (5) 

The percolation point occurs when the non-trivial solution (^J) approaches zero, and this occurs at p c = 1/3. Near p c , 
Td approaches zero linearly, so the transition is second order and the order parameter exponent (3 = 1. 

In order for the problem to lie in the "rigidity percolation" class, there must be at least two degrees of freedom per 
node i.e. g > 2. However when b — 1, if a — z — 1 = 2, then p c = 1, as all bonds must be present to order to transmit 
rigidity. Thus the simplest non-trivial case is g — 2, a — 3 and 6=1, which we now treat. 

Rigidity transition for g — 2, a — 3 and 6=1 
^From the first of Eqs. (2), we have, 

To=p(To 3 + 3T 2 (l-To)) (6) 
Of course there is always the trivial solution To = 0. In addition, Eq. (0) implies 

3±V(9-8/p) 

T ° = 4 W 

To ensure that To = 1 when p = 1, take the positive root. The new feature here is that the square root is negative for 
p < p c = 8/9, so this root is unphysical below p = 8/9. For p < p c , the only remaining real root is Tq = 0, so there is 
a first order jujsip in To at p c = 8/9. The magnitude of this jump ATo = 3/4. Note also that on approach to p c from 
above, we findllH 

Ta-Z/i-ip-pJ 1 / 2 (8) 

which illustrates the singular corrections to the first order jump in To. This interesting behavior seems usual for both 
bootstrap percolation and for rigidity percolation. From the second of Eqs. (2), we have, 

Ti =3pT (l-T ) 2 , (9) 

which has the two solutions, T\ = and the result found by subsituting Eq. (7) for To into Eq. (9). There is thus 
a first order jump in T\ at the same p c as that found for Tq. The size of this jump ATi = 1/8. Note that T\ is zero 
at p = 1, so Ti rises from zero as p is decreases, and peaks at p = p c Since T2 = l — Tq — Ti, all components of the 
vector order parameter are first order, and all of them have a singular correction near p c as a consequence of Eq. (8). 

Order of the transition for general g, a, b = 1 
In the first of Eqs. (2), there is always the trivial solution Tq = 0. After removing that, the following equation holds. 



l=pEu) T o _1 (l-ro) a - fc (10) 

k=g ^ ' 



If g = 1 (connectivity percolation) , there is always a constant term on the RHS of this equation, and this allows a 
real solution for arbitrarily small To, and hence the transition is second-order. However, if g > 2, the constant term 
on the RHS is absent and the equation cannot be satisfied for an arbitrarily small real Tq. Thus there must be a first 
order jump in To for any z > g > 2. It is possible to solve Eq. rtiO) to find p c explicitly in the case g = a — 1, in which 
case the first order jump has magnitude ATo = 1 — l/( a ~ 1) 2 E-j. However in general we resort to numerical methods. 
Before describing the numerical results, we first introduce a matrix method which allows us to treat general g, 6, a 
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B. Site diluted Bethe lattices for arbitrary g,a,b 



It is possible to generalise the Bethe lattices described above to cases where more than one bond connects neigh- 
bouring nodes. In the case of site dilution, removing a site removes all of the b bonds that enter that site from a 
neighbour. In contrast bond dilution removes one bond at a time and must be treated differently(see later in this 
section). Returning to the site dilution case, note that if b > g, rigidity is transmitted across the tree as soon as 
connectivity percolation occurs. This is because any one connection between two nodes with b > g ensures trans- 
mission of rigidity to the newly added node, provided of course that the prior node is also rigid with respect to the 
boundary. Thus if b > g, there are only two possible states for each node: rigidly connected to boundary and not 
rigidly connected to the boundary, and the model is "trivially" in the connectivity percolation class. In contrast, 
if there are fewer than g bars connecting two nodes, more interesting node states are possible, and we must again 
consider the full set To, ....T g , which allow the possibility of partial transmission of rigidity. We now develop a matrix 
method to treat the non-trivial cases 1 < b < g. 

Consider adding a site to a branch of co-ordination a. We label the sites at the previous level i = l,...,a (for 
example the lower ellipse in Fig. lc would have label i = 1). Each of these nodes may have U = 0, 1, g degrees of 
freedom with respect to the border (for example the lower ellipse in Fig. lc has l\ degrees of freedom with respect to 
the border). 

We start by adding a "free body" to the tree, so it has g degrees of freedom with respect to the boundary. However, 
when we add the new higher level body to the tree, we also add ab bonds. But not all of the bonds that are added 
are "useful" in reducing the number of degrees of freedom of the newly added body with respect to the border. For 
example, if a lower level node already has g degrees of freedom with respect to the border, no matter how many bonds 
connect it to the higher level body, it does not produce any constraint of the newly added body with respect to the 
boundary. Therefore we must define the "number of useful bonds" , u, which lies along any sub-branch. If a lower level 
body has zero degrees of freedom with respect to the border, then every bond is "useful" . If the lower level body has 
1 degree of freedom with respect to the border, then the first bond that is added does not constrain the newly added 
node, so that only b — 1 of the bonds are useful. In general if a body has i degrees of freedom, only u — b — i of the 
added bonds are useful in producing constraint in the higher level body. Thus the probability Q u that a sub-branch 
has u useful bonds is given by, (note that since we are considering 1 < b < g, Q g = 0) 

n i Tb - u foru = l,---,6 

Now each sub-branch adds Ui constraints to the newly added body, so the total number of constraints on the newly 
added body is YLt=i u i- Thus the probability that the new node have k degrees of freedom is, 

9 9 9 a 

z 1= o; 2 =o i a =o i=\ 



I, , , : -pEE •- E T^.^Jig - A - f>) 

i 1= 0/ 2 =0 l a =0 i=l 
9-1 

T g = l-J2 T i (12) 

1=0 

Where 9 and 5 are the step function and delta function respectively. 

For numerical purposes, a more convenient way of representing these equations is to add the a sub-branches one 
at a time using a matrix method. We define the vector T L = (T L ,T± , T 2 L , ...,Tg) to denote the probability that the 
newly added body be in one of its possible constraint "states" after the addition of L sub-branches (L = 1, 2... a). If 
we have a free node it has g degrees of freedom so before the addition of any sub-branches, T° = (0,0,0..., 1). We 
then have the recurrence relations, 

T i+1 = T L (T + T 1 + .. + T b ) + T?(T + T 1 + .. + T 6 _i) + .. + T^T a (13) 

and, for I = 1, 2..., g, 

t l+i = T L {Tb + n+i + + Tg) + T L +iTbl + Tl L +2 T b _ 2 + .. + T t L +b T . (14) 
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Eqs.(13) and (14) may be put into matrix form, so that 

f L+1 = MT L = (M) a f° 

with 



(15) 



where, 



/l 01 02 

a T fc _i 
0a 



M 



\o 



T, 



6-1 



0b 
Ti T 



\ 





a Ti 



6-1 

a ) 



and, 



i=b 



b-k 



0k = Y, T i- 



(16) 



(17) 



Finally, we must include the possibility that the site is present or absent, so the probability vector obeys, 

f = p(M) a f" + (l-p)f°. (18) 

As before, the LHS of Eqs. (18) is the probability vector at the next level of the tree in terms of the probabilities at 
the lower levels (which are in the matrix M). 

A little algebra shows that Eqs. (18) reproduce the 6=1 equations (Eqs. (2)) as they must. We illustrate the 
matrix method with a special case (b ^ 1) which is analytically solvable. 

A non-trivial solvable case, a — 2, g = 3,b = 2 
For a, b = 2, g = 3 Eqs. (18) yield, 



To 


+ Tj 


To 





T 2 


+ T 3 


Ti 


To 







T 2 + T 3 


Ti 










To 




i+a-f). 



The first two of these equations yield, 



and 



To = p(T 2 + 2T 1 r ) 



Ti = p(2T X + Tj 2 ). 



(19) 



(20) 



where X = T 2 + T 3 . Since the sum of the T's is one, we have X = 1 — T — Ti and this with Eqs. (19) and (20) yields, 



3T 2 - 4(2 - l/p)T + 1/p 2 = 0, 



(21) 



Solving for T yields, 



G 



T = 



(4p-2) + 2V((2p-l) 2 - 
3 P 



3/4) 



(22) 



Then the argument of the square root becomes negative for p < p c given by, p c = (1 + \/3/2)/2 ~ 0.933, so that 
AT = 0.619. 



Numerical results for general b, g, a 

Results of iterating the matrix Eqs. (18) are presented in Figs. 2-4. Fig. 2a illustrates that for g < b, the problem 
reduces to the connectivity percolation case. The transition is second order and only two components of the vector T 
(To and T g ) are finite. In contrast, when ba > g > b (see Fig. 2b), all of the components of T can be finite, although 
all of them are singular at the same percolation point. This figure also illustrates the point that the rigidity transition 
is first order and we have selected this case to illustrate the fact that in some case the rigidity transition is weakly 
first order. 

In Fig. 3, we illustrate the dependence of rigidity percolation on the co-ordination number a. In the case we choose 
here, g = 2, b = 1, the transition is always strongly first order. The behavior near p = 1 is typical of site dilution on 
any lattice, because the leading term in the probability that a site is not rigid with respect to the boundary, is just 
the probability that the site is absent, i.e. 1 — p. As a increases, the point at which T breaks away from 1 — p tends 
to p = as intuitively expected. 

If we start from a rigid border, it is easy to verify that the transmission of rigidity depends on a and the ratio b/g. 
in the limit b/g — > 1, we have rigidity percolation, while if b/g — > 0, the transition is at p = 1 and is completely first 
order. Using trees, we are able to probe various values of b/g and we present results for p c (a, b/g) in Fig. 4. It is seen 
that for all cases, p c ~ ^1 for a — > oo. We also find that for any b/g < 1, the transition is first order, and the size 
of the first order jump increases smoothly as b/g decreases. 

^From the site dilution problem, we conclude that the rigidity transition is always first order, except in cases where 
it trivially reduces to connectivity percolation. However, there appears to be a square root singularity superimposed 
on the first order jump in T . However, on site diluted lattices with b < g, the only rigid clusters are those which 
are attached to the rigid border. In contrast in bond percolation it is possible to have internal rigid clusters, and the 
cases b > g are non-trivial. Thus we now describe calculations for the transmission of rigidity in bond diluted trees. 



As for the site diluted case, we define the vector T = (T , Ti, T 2 , T g ). Now, if there is a total of b possible bonds 
between two nodes, and if each is present with probability p, then the probability q(k) that k bonds are actually 
present is 



Since the nodes have g degrees of freedom, at most g independent bonds can connect two nodes. If k > g bonds 
connect two nodes, k — g of them will be redundant and the two nodes will form part of a cluster that is internally 
rigid. Any number of bonds in excess of g does not add to the number of independent constraints. Therefore the 
probability qk that k independent bonds are present between two nodes is, for general g and b, 



As in the site dilution case, these k bonds are not all "useful" in transmitting constraint from the boundary unless 
the sub-branch along which they lie is at least partially constrained. In fact if the lower level node has i degrees of 
freedom with respect to the boundary, only k — i of the bonds connecting that node to the newly added node actually 
impose contraint. Clearly if k < i, the branch imposes no constraint (with respect to the boundary) on the newly 
added node. We thus define the useful bonds u = k — i , because they are able to propagate contraint outward from 
the boundary. The probability Q u for a branch to have u useful bonds on it is then given by, 



C. Bond-diluted Bethe lattices (general g, a and b) 




(23) 




for k < g 
for k = g 
for k > g 



(24) 




for u = 1, • • • 
for u = 



,9 



(25) 



Now taking a such sub-branches, the total number U of useful bars is 
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U = J2 u k (26) 



fe=i 



If C/ > <7, then the new node body will be rigid. Otherwise it will have k = g — U degrees of freedom. Formally we 
then write 



g g 



T f=J2J2'-- 2 Q U iQ« a ---Q«,_ 1 *(/,«i,«2,-",««-i) (27) 

K1-O112-O u z -i— 



where 



f <5((7 - (.g - /)) for < / < g 

3>{f,g,z,U) = \<m& (28) 
I 0(tf - ff) for / = 

Where as in the site case, we have used the step function and the kronecker delta to ensure that the constraint 
counting is correct. 

As for the site diluted case, we can write the Eqs. (27) in matrix form. 

f L+1 = MT L (29) 

with 



M 



/I (Ql+Q2 + --- + Pg) {Q2 + Q 3 + --- + Qg) ■ 
Qo Ql 

Qo 



Vo 











Q 9 \ 

Qg-l 

Qg-2 

Qo J 



Again starting from a bare node with T° = {0, 0, • • • 0, 1}, and after connecting a legs, we get the desired probabilities 

as 



f = {M) a f° 

To illustrate the matrix method for the bond case, we again do a solvable case with b 7^ 1. 
A non-trivial solvable case a = 2, g = 3, b = 2 
,iFrom Eqs. (23) and (24), we have, 

(50,91,52,93) = ((1 -£>) 2 ,2p(l -p),p 2 ,0). 

Then from Eq. (25), we have, 

(Qo, Qi, Q2, Qs) = (1 - (P 2 + 2p(l - p))T - P 2 T U 2p(l - p)T +p 2 r!,p 2 T , 0). 
Using these expressions in the matrix equation (29), we have, 



(30) 



(31) 
(32) 



/T \ (l (2p-p 2 )T +p 2 T 1 p 2 T Q \ /°\ 

Ti I 1 - (2p - p 2 )T - p 2 T! 2p(l - p)T + p 2 Tx p 2 T 

T 2 l-(2p-p 2 )To-p 2 T! 2p(l - p)T + p 2 Tx 

Vt 3 / VO l-(2p-p 2 )T -p 2 r 1 / \lj 

^From the first of these equations, we find 

To = P 3 T {(4 - 3p)T + 2pT!}, (33) 

while the second implies 
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Ti = 2p 2 T (l - (2p - p 2 )T - p 2 ^) + (2p(l - jj)T + p 2 Tx) 2 . (34) 
Solving Eqs. (33) and (34) for To gives the trivial solution To = 0, and, 

rp _ (8p 3 - 12 P + 8) + V((8p 3 - Up + 8) 2 - 12p 2 ) 

To 6^ ( 35 ) 

This again becomes imaginary at the rigidity threshold, which we find to be p c = 0.918, and the first order jump in 
T is, AT = 0.629. 

Numerical results for general 6, g, a 

First we note that the for 6 = 1, site dilution and bond dilution are the same, provided we make the transformation 

Psite * Pbond 

and T s ite — Pbon dTbond, thus we focus attention on 6 > 2. 
We present numerical results for bond diluted trees in Figs. 5 and 6. In Figure 5, we show that even when b >> g 
and many internal rigid clusters can exist on the trees, the rigidity transition remains first order. In fact, we have not 
found any values of g or 6 for which the bond diluted trees are second order, except the trivial case g = 1. However 
the rigidity transition is weakly first order for b/g — > oo. A second interesting feature of Fig. 5 is the non-monotic 
behavior of T±. Nevertheless on all of the trees we studied, the rigidity transition is unique and first order. As in Eq. 
(8), there appears to be a singular behavior superimposed on the first order jump in Tq. On the bond diluted trees, 
the percolation threshold depends on all three parameters g,b,a, nevertheless there is a simple behavior in the large 
a limit (see Fig. 6), so that p c ~ G(g, b)/a for a — > oo. 



IV. MECHANISM AND COMPARISON WITH OTHER THEORIES 



A mechanism for first order rigidity 

The mechanism for the first order rigidity transition is illustrated in Fig. 7a for an a = 2, g = 2, b = 1 tree and in 
Fig. 7b for the bond-diluted triangular lattice. In these figures, we have presented a rigid cluster and have indicated a 
bond which we then remove. On removal of the arrowed bond, both of the rigid clusters "break" up into more than 2 
rigid subclusters. In Figure 7a, removal of the arrowed bond leads to 6 rigid subclusters, while in Fig. 7b, removal of 
the arrowed bond leads to the formation of 4 rigid subclusters. In both cases wc arc referring to clusters of mutually 
rigid bonds. In contrast, in connectivity percolation, removal of a "cutting" or red bond leads to the break-up of the 
system into two subclusters. On large rigid clusters, the removal of a "cutting" or red bond usually leads to formation 
of many subclusters, and this "cluster collapse" provides a mechanism for a first order rigidity transition. However 
it does not ensure a first order transition, as it depends on how many clusters are formed when a cutting bond is 
removed. In reverse the phenomenon of cluster collapse is "cluster-freezing" in which there is a sudden jump in the 
average cluster size as many clusters suddenly become mutually rigid (For example by replacing the arrowed bonds 
in Fig. 8). It is likely that these ideas can be used to develop scaling arguments for the amount of cluster-collapse 
required for there to be a first order rigidity transition, and we are currently working in that direction. 

Comparison with constraint counting methods 

For simplicitiy, consider first bond percolation for which the argument is simplest. On a regular lattice, there are 
N nodes of co-ordination z, with each node having g degrees of freedom and with 6 bars connecting each pair of 
nodes. Now dilute the bars of the network, with p the probability rtjaat any one bar is present. Then "on average", 
the number of degrees of freedom, fN, that remain at dilution p isc3, 

fN = Ng - pbzN/2 + B, (36) 

where the factor of 1/2 is due to the fact that each bar is shared between two nodes. B is the number of bonds that 
are "redundant" in that they are in regions of the lattice which would be rigid even if they were removed. The mean 
field approximation reduces to assuming B = 0, so that / = g — pbz/2 and thus / approaches zero at p c = 2g/bz. 
This counting procedure is slightly modified on trees, as the border is rigid so every bond which is next to but lower 
than a node in the tree contributes to the rigidity of that node (the bonds are not "shared" as on a regular lattice). 
In this case, the constraint counting is 

fN = Ng - pbaN + B. (37) 

Thus we have the same expression as in Eq. (36), with the replacement a(tree) ~ z/2(regular lattice). If we again 
assume that B — 0, we find, p c (B = 0) = g/(ba). This estimate is grossly in error when compared with the actual 







results for trees (see Fig. 6). Clearly the stronger the first order transition, the more in error the constraint counting 
mean field theory becomes. 

Global constraint counting I— . 

It has been observed that in two dimcnsionsliJ, although the number of floppy modes is always continuous, the 
second derivative of that quantity is singular. This is based on counting the number of degrees of freedom in the 
whole lattice. If we do a similar calculation on trees, the surface bonds dominate, nevertheless it is interesting to see 
what the results are. Thus we have done a calculation which keeps track of the number of redundant bonds on the 
trees for all levels going outwards from a rigid boundary. We have done the calculation for bond diluted lattices with 
6=1. In that case, the number of redundant bonds I levels away from the boundary is given by, 

B, = a L ~ l f> - g) (fj (pT'-y (1 - pT^ 1 )""* (38) 

k=g ^ ' 

L is the total number of levels in the tree. The total number of redundant bonds in the tree is, 

L 

B = ]>>,. (39) 
i=i 

^From global constraint counting, we then have, 

f = g-pa + B/N s (40) 

N s = a L I (a — 1) is the number of sites on the L level tree. Results for /, df /dp and d 2 fdp 2 are presented in Fig. 8. 
It is clear from these calculations that there is no singular behavior in the second derivative of / on trees. However, 
there is a peak in the second derivative, but at a value of p considerabley less than p c . 



V. CONCLUSIONS 



We have shown that it is straightforward to develop and analyse tree models for the transmission of rigidity from a 
rigid border. In order to analyse these models we must, in general, consider the transmission of "partial" rigidity, as 
partially rigid structures may lead to rigidity higher up the tree. Some of the main conclusions of our calculations are 

1. Except for some "trivial" cases which are equivalent to connectivity percolation, the rigidity transition in these 
systems is first order. However there may be a singular piece superimposed upon the first order transition in the 
infinite cluster probability, as was explicitly demonstrated in some special cases (see e.g. Eq. (8)). 

2. Constraint counting mean field theory which ignores redundant bonds is qualitatively incorrect for trees. This 
method does not describe correctly the nature of the rigidity transition. It can also grossly underestimate p c , especially 
if the transition is strongly first order. 

3. We have defined a vector order parameter which describes the number of degrees of freedom two points have 
with respect to each other. Although there is the possibility of multiple phase transitions with such a vector order 
parameter, we find that there is only one transition on trees. 

4. The number of floppy modes and its first and second derivatives are non-singular, probably due to the dominance 
of surface bonds on trees. 

5. Bootstrap percolation and rigidity percolation are exactly the same on 6 = 1 trees, but different on regulan-iattices. 
It is not clear, at least to these authors, to which case (if either), the current continuum field theory appliestil 

Taken together with new numerical results in two and three dimensions@lll@, there is now quite strong evidence 
that the rigidity transition on random lattices is often first order, in contrast to the large number of earlier papers 
which have assumed the opposite. However it is important to emphasize that the new work using exact constraint 
counting is correct for random lattices while the earlier work was for regular lattices. It is still an open question as to 
whether these two cases are qualitatively different. 
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1 The geometry of trees, a) A z — 5, b — 1 tree; b) One branch of the tree of a); c) One branch of a b = 3 tree. 

2 Rigidity percolation of site diluted trees: a) a — 4, g = 3 and 6 = 3. The infinite cluster probability and the 
probability T 3 are plotted. In this case the behavior is the same as connectivity percolation, so p c — l/a and the 
transition is second order, with /3 = 1; b) a = 5, g = 3 and b = 2. T , T\, T 2 and T 3 are plotted. All are first order 
and singular at the same p c 



3 Tq for g = 2 and a = 1 for various a. The transition is always first order 

4 p c as a function of b/g and a. From the top, the curves are for b/g = 1/6; 1/3; 1/2; 2/3; 1 



5 Rigidity percolation for a bond-diluted tree with a = 2, g = 2 and b = 40. The transition is close to second order 
and there is an interesting non-monotonic behavior in T\ 

6 p c for bond-diluted trees. Curves are for (from the top) g — 3, b = 1; g = 2,6 = 1; q = 6,b = 5; g = Q, 6 = 10; 
.9 = 2,6=10 

7 The effect of removing a bond on the cluster size distribution, a) Removing the arrowed bond from this rigid 
cluster leads to 6 separate rigid clusters, b) Removing the arrowed bond from this connected cluster leads to 4 separate 
rigid clusters. 

8 Floppy modes on a bond-diluted tree with a — 6, g = 3 and 6 = 1. The number of floppy modes per site is 
continuous as are its first (/' and second /"derivatives. 
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